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We present a determination of the CKM parameter |Kjs| based on new measurements of the 
six largest Kl branching fractions and semileptonic form factors by the KTeV (E832) experiment 
at Fermilab. We find \Vus\ = 0.2252 ± O.OOOSktoV ± 0.0021oxt, where the errors are from KTeV 
measurements and from external sources. We also use the measured branching fractions to determine 
the CP violation parameter [j7+_| = (2.228 ± 0.005ktcV ± O.OOQoxt) x 10"^. 



PACS numbers: 12.15.Hh, 13.25.Es, 13.20.Eb 

The Gabibbo-Kobayashi-Maskawa (GKM) matrix |l|, 
1^ describes the charged current couplings of the u, c, and 
t quarks to the d, s, and b quarks. The first row of this 
matrix provides the most stringent test of the unitarity of 
the matrix. Gurrent measurements [2| deviate from uni- 
tarity at the 2.2 sigma level: l-(|KdP + |KsP + |Kb|^) = 
0.0043 ± 0.0019. \Vus\, which contributes an uncertainty 
of 0.0010 to this unitarity test, has been determined from 
charged and neutral kaon semileptonic decay rates. This 
determination is based on the partial width for semilep- 
tonic K decay, Tkis'- 



In this Letter, we present a determination of \Vus\ by 
the KTeV (E832) experiment at Fermilab based on mea- 
surements of the Kl — > ir^e^i' and K^ — > ir^ jM^v par- 
tial widths and form factors. These measurements are 
described in detail elsewhere |j, |a| ; a brief summary is 
given here. Our |Kis| determination also makes use of a 
new treatment of radiative corrections 6] . 

To determine the Kl 
partial widths, we measure the following five ratios: 



n e^v and Kl — > it jX^v 
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where i refers to either e or ^, Gp is the Fermi constant, 
Mk is the kaon mass. Sew is the short-distance radia- 
tive correction, 5^^ is the mode-dependent long-distance 
radiative correction, /+(0) is the calculated form fac- 
tor at zero momentum transfer for the (.v system, and 
/^ is the phase-space integral, which depends on mea- 
sured semileptonic form factors. C^ is 1 (1/2) for neutral 
(charged) kaon decays. The current PDG determination 
of \Vus\ is based only on K -^ -nev decays; K — > -niiv 
decays have not been used because of large uncertainties 
in/j^. 
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where internal bremsstrahlung contributions are included 
for all decay modes with charged particles. Since the six 
decay modes listed above account for more than 99.9% of 
the total decay rate, the five partial width ratios may be 
converted into measurements of the branching fractions 
for the six decay modes. The Kl lifetime is then used 
to convert these branching fractions into partial widths. 
The branching fraction measurements also can be used to 
determine the GP violation parameter \r]^ p = ^{Kl -^ 
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FIG. 1: The KTeV (E832) detector. 



The semileptonic form factors describe the distribu- 
tion of t, the momentum transfer to the Iv system. 
This t dependence increases the decay phase space in- 
tegrals, /^ and I^, by about 10%. We use the foUow- 
ing parametrization for the two independent semileptonic 
form factors: 



/+W = /+(o)(i + a;^ + U" *' 
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where /+ (0) is obtained from theory, and we measure A^ , 
A" , and Aq. 

In principle, the form factors can be measured directly 
from the t distribution. The undetected neutrino and 
unknown kaon momentum, however, result in a twofold 
ambiguity in the reconstructed value of t. To avoid sys- 
tematic uncertainties associated with this ambiguity, we 
use a technique based only on components of particle mo- 
menta measured transverse to the kaon momentum J5|. 

The KTeV experiment (Fig. ^ and associated event 
reconstruction techniques have been described in detail 
elsewhere |3|. An 800 GeV/c proton beam striking a 
BeO target is used to produce two almost parallel neu- 
tral beams. The regenerator beam, which includes Ks, 
is not used in this analysis; the vacuum beam provides 
K]^ decays used for these measurements. A large vac- 
uum decay region surrounded by photon veto detectors 
extends to 159 m from the primary target. Following the 
vacuum region is a drift chamber spectrometer, trigger 
hodoscope, pure Csl electromagnetic calorimeter, and a 
muon system consisting of scintillator hodoscopes behind 
4 m and 5 m of steel. The analyses presented in this Let- 
ter make use of the detector calibration and Monte Carlo 
simulation from the KTeV e'/e analysis [3. 

Simple event reconstruction and selection may be used 
to distinguish different kaon decay modes from each 
other, and to reduce background to a negligible level for 
all decay modes. The reconstruction of charged decay 
modes {Kl — *■ n^e^v, K^ —^ tt jjl^v^ Kl — *■ tt'^tt^tt^ , 
and Kl -^ tt+tt^) begins with the identification of two 



oppositely charged tracks coming from a single vertex. 
Note that for the K^ -^ 7r"'"7r~7r° decay, we choose not to 
reconstruct the 7r° -^ 77 decay to reduce the acceptance 

uncertainty in the r_| o/^Ke3 ratio. 

The charged decay modes are separated from each 
other on the basis of particle identification and kinematic 
requirements. To select Kl -^ ir^e^v decays, the elec- 
tron is identified using the calorimeter energy measure- 
ment (E), combined with the spectrometer momentum 
(p). Kl — *• 7r+7r~ is separated from other two-track de- 
cays based on the two-track invariant mass, m^7r, and 
the square of the two-track momentum transverse to the 
Kl direction, p^. To isolate K^ — > 7r+7r~7r° we use an 

additional variable, fc_| Qj described in [j|. Note that 

the Kl ~^ ■k'^^j^v background to each of these decay 
modes is suppressed by rejecting events with hits in the 
muon system. Two-track events that are not identi- 
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fied as Kl -^ tt e 
are selected as Kl -^ ir'^ii^y candidates. To reduce the 
acceptance uncertainty in the Tkiis/^KcS ratio, we do 
not require a signal in the muon hodoscope to identify 
Kl — > -K^p^v decays. 

The reconstruction of the Kl -^ tt^tt*' and 
Kl — > Tr'^Tr'^Tr'^ decay modes, where tt^ -^ 77, is 
based on energies and positions of photons measured in 
the Csl electromagnetic calorimeter as described in \(\. 
Exactly four (six) clusters, each with a transverse profile 
consistent with a photon, are required for Kl — > tt^tt" 
{Kl — > 7r°7r°7r°). Photons are paired to reconstruct two 
or three neutral pions consistent with a single decay 
vertex. 

All reconstructed decay modes are required to have 
kaon energy, Ek^ between 40 and 120 GeV, and decay 
position, z, between 123 and 158 m from the target. For 
the reconstruction of semileptonic and Kl — > tt'^tt~tt^ 
decays, there is a missing particle {i^ or tt"), which re- 
sults in multiple kaon energy solutions. Each of these 
solutions is required to be in the accepted range. 

After all event selection requirements and background 
subtraction, we have between 10^ and 10^ events per de- 
cay mode. The background is 0.7% for Kl -^ tt^tt" and 
much smaller for the other decay modes. After correct- 
ing each of the ratios for acceptance differences between 
numerator and denominator, we find the partial width ra- 
tios given in TablcQ] The precision is 1.2% for Tooo/^KeS 
and about 0.5% for other ratios. A comparison of data 
and MC z-vertex distributions for the semileptonic de- 
cay modes (Fig. El demonstrates the quality of the MC 
simulation used for the acceptance correction. 

The five partial width ratios may be combined to de- 
termine the branching fractions shown in Table |n] S . 
Using the PDG average for the neutral kaon lifetime |3], 
tl = (5.15 ± 0.04) X 10~^ s, our branching fraction mea- 
surements correspond to the partial decay widths shown 
in Table ni 

Figure |3| shows a comparison of the KTeV and PDG 



TABLE I: Measured partial width ratios. The first error 
is statistical and the second systematic. The five statistical 
errors are independent; correlations among the systematic er- 
rors are treated in 41. 



Decay Modes 



rooo/rife3 

r+-o/rKe3 

r+-/rxe3 

roo/rooo 



Partial Width Ratio 



0.6640 ± 0.0014 ± 0.0022 

0.4782 ± 0.0014 ± 0.0053 

0.3078 ± 0.0005 ± 0.0017 

(4.856 ± 0.017 ± 0.023) x 10"^ 

(4.446 ± 0.016 ± 0.019) x 10"^ 



TABLE IL Kl branching fi-actions and partial widths iVi). 
Correlations among uncertainties in these measurements are 
given in |4|. 



Decay Mode 



Branching fraction 
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865 ±0.010) X 10"^ 



0.7897 ± 0.0065 
0.5244 ± 0.0044 
0.2431 ± 0.0023 
0.3777 ± 0.0045 
(3.835 ± 0.038) x 10"^ 
(1.679 ±0.024) X 10"^ 
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FIG. 2: (a) Comparison of the vacuum beam z distributions 
for data (dots) and MC (histogram); (b) Data-to-MC ratios 
with a linear fit. 



values for the six branching fractions. The new KTeV 
measurements are on average a factor of two more pre- 
cise than the current world average values, but are not in 
good agreement with these averages. Compared to the 
PDG fit 0, the KTeV measurement of B{Kl -^ n^e^v) 



is higher by 5%, B{Kl -^ 7r°7r°7r°) is lower by 8%, 
B{Kl -^ vr+TT-) is lower by 5%, and B{Kl -^ Tr^Tr") is 
lower by 8%. Our measurements of B{Kl — > 7r='=^±i') 
and B{Kl -^ 7r+7r~7r°) are consistent with the PDG fit. 
A detailed comparison between the KTeV measurements 
and previous results is given in |j]. 



Using the measured branching fractions for 
Kl -^ vr+vr^ and Kl — ^ tt'^tt" together with 
tl, ts, Re{e'/e), and B{Ks -^ n£iy), we de- 
termine the CP violation parameter \ri^ | = 

(2.228 ± 0.005ktoV ± 0.009ext) x IQ-^; most of the 
external error results from the uncertainty in tl- Our re- 
sult is 2.6% lower than the PDG average. A comparison 
of |?7_| I determinations is given in J4f|. 

The /+ (t) form factor is measured in both semileptonic 
decay modes; the effect of fo{t) is proportional to the lep- 
ton mass, so it is only measured in Kl — ^ tt fi^iy decays. 
The measured parameters for the semileptonic form fac- 
tors are A^ = (20.64±1.75) x IQ-^, A^^ = (3.20±0.69) x 
10-3, and Ao = (13.72±1.31) x IQ-^. The corresponding 
phase space integrals are /^ = 0.15350 ± 0.00105 and 
/^ = 0.10165 ±0.00080, where the quoted errors include 
an additional uncertainty related to the form factor pa- 
rameterization I5j|. Compared to phase space integrals 
based on PDG form factors, KTeV's /J^ and I'^ integrals 
are 1.7% and 4.2% lower, respectively. If we fit our data 
without the A" term, our /^ and I'^ integrals are in- 
creased by 1%, and are consistent with PDG averages 
that use only linear terms. 

To check the consistency of our branching fraction and 
form factor measurements with lepton universality, we 
compare Gp for the two decay modes by taking the ratio 
of Eq. ^for Kl — > n^fi^i^ and Kl -^ -K^e^v: 
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FIG. 3: Kl branching fractions measured by KTeV (dots) 
and from PDG fit (open circles). 
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Many common uncertainties cancel in this ratio. The 
ratio of radiative corrections is calculated to be (1 ± 
SK)/i^ + ^k) = 1-0058 ± 0.0010 ?], the ratio of the 
phase space integrals is I'^/I'x = 0.6622 ± 0.0018, and 
^Kiiz/^KeZ is from Table Ull The resulting ratio of cou- 
plings squared is {G%/G%f = 0.9969 ± 0.0048, con- 
sistent with lepton universality. The same ratio calcu- 
lated from PDG widths and form factors is (Gp/Gp)^ — 
1.0270 ± 0.0182. Note that the 0.5% uncertainty in our 
universality test is much smaller than the 5% difference 
between the KTeV and PDG values of TKfis/^KeS- 

The measured partial widths and phase space inte- 
grals for semileptonic decays can be combined with theo- 
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FIG. 4: Comparison of tlie KTeV measurement of |Kjs|,f+(0) 
with Broolcliaven E865 |lOll . PDG, and also with determina- 
tions oi f+{0){l-\Vudf -\Vub^^^ based on different theoret- 
ical calculations of /+(0) p. Illl [l2 . U ^ . K^ measurements 
have been divided by 1.022 ± 0.005, the ratio of /+(0) for 
charged and neutral kaons 1141 . For Kl measurements, the 
uncertainties are mainly from tl- PDG refers to our eval- 
uation based on PDG partial widths and form factors. For 
/+(0)(1 — IKjdP ~ IKjfcl^)"'^ , the inner error bars are from 
/+(0) uncertainty; the total uncertainties include the \Vud\ 
and \Vub\ errors. 



retical corrections to calculate | Vus \ using Eq. Q] The 
short-distance radiative correction. Sew = 1.022 [sj, 
is evaluated with a cutoff at the proton mass. The 
long-distance radiative corrections are taken from [g: 
(5^ = 0.013 ± 0.003 and <5^ = 0.019 ± 0.003. For /+(0), 
we use the same value used in the PDG evaluation of 
\Vus\- /+(0) = 0.961 ±0.008 [3. 

The resuhing values of |Ks | are 0.2253±0.0023 for K^s 
and 0.2250 ± 0.0023 for K^^, where the errors include an 
external uncertainty of 0.0021 from /+(0), the Kl hfe- 
time, and radiative corrections. Assuming lepton univer- 
sality, we average the K^ 
results (accounting for correlations): 



TT e^i' and Kl -^ tt ji^v 



\Vus\ = 0.2252 ± O.OOOSktcV ± 0.0021e 



(9) 



The KTeV error comes from uncertainties in the KTeV 
branching fraction and form factor measurements. 

To compare our result with previous charged and neu- 
tral kaon measurements, we use the product of \Vus\ and 
/+(0) rather than \Vus\ to avoid significant common un- 
certainties from /_|_ (0) . Figure ^ shows a comparison of 
our measurement of 



|Ks|/+(0) = 0.2165 ±0.0012 



(10) 



with values from the PDG and Brookhaven E865 |lO| . 
Our value of |Kis|/-i-(0) is inconsistent with previous 
Kl determinations, but is consistent with K^ results 
(both E865 and earlier measurements). The figure also 



shows /+(0)(1 - IKdP - IKbP)^^^, the expectation for 
/+(0)|Kis| assuming unitarity, based on \Vud\ = 0.9734 ± 
0.0008, iKbl = (3.6 ±0.7) x 10^^^ and several recent cal- 
culations of /+(0). Our value of \Vus\ (Eq. O, based 
the Leutwyler and Roos calculation of /+(0), is con- 
sistent with unitarity: 1 — (|Kid|^ + |KisP + iKifoP) = 
0.0018 ±0.0019. For other calculations of /+(0), the con- 
sistency with unitarity ranges from 1 to 1.7 sigma, as 
shown in Fig. 0] Our improved form factor measure- 
ments ma y h elp to reduce theoretical uncertainties in 

./+(o) mEHIa. 

In summary, KTeV has made improved measurements 
of the six largest Kl branching fractions and the semilep- 
tonic form factors. We use these results to deter- 
mine |?7+_| = (2.228 ± 0.010) X 10^^ ^nd |Kc|/+(0) = 
0.2165 ±0.0012. Using /+(0) = 0.961 ±0.008 0, we find 
Kisl = 0.2252 ± 0.0022, consistent with unitarity of the 
CKM matrix. 
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